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At zero temperature, two-site dynamical mean field theory is applied to the Dynamic Hubbard 
model. The Dynamic Hubbard model describes the orbital relaxation that occurs when two electrons 
occupy the same site, by using a two-level boson field at each site. At finite boson frequency, the 
appearance of a Mott gap is found to be enhanced even though it shows a metallic phase with 
the same bare on-site interaction U in the conventional Hubbard model. The lack of electron- 
hole symmetry is highlighted through the quasi-particle weight and the single particle density of 
states at different fillings, which qualitatively differentiates the dynamic Hubbard model from other 
conventional Hubbard-like models. 



PACS numbers: 

I. INTRODUCTION 

Dynamic Hubbard models^ represent a new class of 
model Hamiltonians that describe the modification of the 
electronic wavefunction that occurs when an atomic or- 
bital becomes occupied by more than one electron. The 
key difference between dynamic and standard Hubbard- 
like models can be understood by considering a Helium 
atom. Helium has two electrons in the Is shell and the 
strength of the interaction between them is comparable 
to the electron-ion interaction. As a consequence, the 
two-electron wave function cannot be simply represented 
by a single Slater determinant formed by the electronic 
wavefunction of the singly occupied orbital (He + ) for 
each electron, since the wavefunction of one electron is 
modified by the presence of the other electron. The ef- 
fect is more pronounced in the negative ion H~ because 
of the weaker attraction between the electrons and the 
nucleus. Conventional tight binding models like the Hub- 
bard model completely neglect this effect since they as- 
sume that the wavefunction for two electrons on a site is a 
simple product of the wavefunctions for a singly occupied 
sitei. 

The Dynamic Hubbard model tries to include the 
physics of orbital relaxation by modulating the on-site in- 
teraction term (Hubbard U) with an auxiliary boson (or 
spin) degree of freedom^ An effective low-energy Hamil- 
tonian for this model modifies the hopping term so that 
the hopping amplitude depends on the electronic occu- 
pation of the sites involved in the hopping process; this 
is known as the correlated hopping model 2 . In this pa- 
per, we focus on the dynamic Hubbard model with an 
auxiliary spin- 1/2 degree of freedom; the two states of 
this pseudospin can be viewed as describing the modi- 
fication of the electronic wavefunction upon double oc- 
cupancy. This kind of model was first suggested two 



decades ago4 The model was studied further through 
world line quantum Monte Carlo (QMC) methods^, ex- 
act diagonalization (ED)£, and an approximate pertur- 
bative analysis starting from a generalized Lang-Firsov 
transformation^ The exact treatments suffer from small 
size effects, but more recently, Bouadim et al.— also stud- 
ied this model using determinant QMC on somewhat 
larger clusters (N = 6x6); however, due to the fermion 
sign problem, they encountered difficulties, especially at 
low temperatures. 

In recent years it has become clear that dynamical 
mean field theory (DMFT) is a valuable way to treat 
the local aspects of both quasiparticles and incoherent 
high energy excitations on the same footing in strongly 
correlated electron systems. 8 DMFT is practical because, 
instead of considering a large lattice model whose Hilbert 
space is exponentially large, one needs to solve merely 
a single impurity model; for this problem various algo- 
rithms for an exact solution already exists Even though 
the single impurity model can be treated numerically 
by QMC, it is still a computationally expensive prob- 
lem. An approximate but effective alternative was pro- 
posed by PotthoffjiS who proposed the so-called "two-site 
DMFT"; this method uses two sites, one for the impu- 
rity and one for the 'bath' of conduction electrons, which 
is readily solved exactly. Using only one bath site ren- 
ders the mapping to the single impurity Anderson model 
(SIAM) approximate, but the self-consistency is now eas- 
ily controlled through the bath parameters. In Ref. [13 
the validity of the two-site DMFT simplification was es- 
tablished for the Hubbard model; for example, it predicts 
qualitatively correctly the existence of a Mott transition 
critical point, and even though it is not a conserving the- 
ory in the sense of Baym and Kadanoff, 11 the violation 
of the Luttinger sum rule, for example, is fairly small. 
Given our interest in establishing qualitative trends for 
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the Dynamic Hubbard model, and the existence of QMC 
results (albeit for small lattices) as a benchmark, we will 
adopt the two-site DMFT approximation to study the 
dynamic Hubbard model. 

There have been some studies of correlated hopping 
models using DMFT, in particular applied to the Falicov- 
Kimball model 12 ' 13 . Because the interaction term in- 
volves two sites rather than one, the self-energy is non- 
local and the formalism becomes considerably more com- 
plicated. Instead with the dynamic Hubbard model con- 
sidered here we can describe the physics of correlated 
hopping in a much simpler way with a single-site self- 
energy, by simply considering the model in the limiting 
case where the interaction becomes non-retarded (large 
wo in the Hamiltonian Eq. (1)). 

As we illustrate below, the two-site DMFT treatment 
of the Dynamic Hubbard model gives semi-quantitative 
agreement with the QMC results. For half-filling and be- 
low, the properties of the Dynamic Hubbard model mimic 
those of the Hubbard model. For example a Mott insu- 
lating phase appears for strong enough on-site interac- 
tion, and, with an attractive on-site interaction, pairing 
is enhanced. In addition, however, electron-hole asym- 
metry naturally arises; this is evident, for example, in 
the dependence of the quasi-particle weight on the elec- 
tron/hole number density. Thus this model captures the 
essential physics that a few electrons in a nearly empty 
electronic energy band can behave very differently from 
a few holes in a nearly filled band. 

The paper is organized as follows: The next section 
will briefly describe the Dynamic Hubbard model with 
an auxiliary spin 1/2 degree of freedom and will pro- 
vide a synopsis of the two-site DMFT approximation. In 
section III we present some numerical results and discuss 
some of the characteristic properties of the dynamic Hub- 
bard model, especially those that differentiate it from the 
simple Hubbard model. In addition, we show some com- 
parisons with the QMC results^ which indicate that the 
two-site approximation works very well. The last section 
IV will summarize our results and suggest directions for 
further study. 



II. MODEL AND METHOD 

We consider here the Dynamic Hubbard Hamiltonian 1 
with a spin- 1/2 degree of freedom in the electron repre- 
sentation: 

-HdHM = tij { c \ a c ja + c }a C iv) ~ M ^ n ia 

<i,j>cr i.a 

+ ^2(uo<7? + guoo-f) + ^(U - 2guj Q a z )n tf n zl . (1) 

i i 

The first term is the electron hopping term; c\ a (ci a ) 
is an electron creation (annihilation) operator at site i 
with spin a. Following Potthoff 1 - we use a Bethe lattice 
with infinite connectivity with nearest neighbour hopping 
only, so that = — t < for nearest neighbours only. 
The parameter t = t* / \fq, with q the connectivity, and 
t* = 1 sets the energy scale. The second term is the 
usual chemical potential term which determines the elec- 
tron filling. The auxiliary spin degree of freedom is given 
in the third term; the two levels have a spacing given 
by loq. The fourth term describes interactions between 
two electrons. In addition to the onsite Hubbard U term, 
there is an additional coupling to the auxiliary spin de- 
gree of freedom. As explained in Ref. [l| and reviewed 
in the next section, this term varies the actual on-site 
repulsion, dependent on the state of the auxiliary degree 
of freedom. 

Dynamical Mean Field Theory has been widely stud- 
ied in a number of correlated fermion systems, and this 
approach has been quite successful, as reviewed in Ref. 
8. In particular, in models where the local dynamics is 
expected to play the most important role (as opposed 
to spatial correlations), the DMFT without the use of 
cluster mcthod o 14 i 15 should be accurate. The Dynamic 
Hubbard model should be ideally suited for these con- 
ditions. DMFT maps a lattice model onto a quantum 
single impurity model through self-consistent conditions; 
in this paper, we consider particularly the single impurity 
Anderson model (SIAM): 
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a a,k=2 a,k=2 
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where d a , dj. are the impurity operators with spin a 
which only act on the single impurity orbital with en- 
ergy ed- This impurity site is hybridized to a bath with 
n s — 1 degrees of freedom through the coupling parame- 
ter 14; the Coulomb interaction (U — 2gojQ(j z ) only occurs 
when two electrons are on the impurity site. The on-site 



Green's function for the lattice model is given by 

p(x) 



G(w) 



dx- 



U! + jJL — S(w) — X 



(3) 



where S(to) is the local self-energy in infinite dimensions 
(q — > oo) and p{x) is the free density of electron states 
for a Bethe lattice: 
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p(x) 



2irt* 



(4) 



Note that a 'momentum' dependent spectral Green func- 
tion is implied by Eq. ([3]), and, for momentum corre- 
sponding to bare energy x, is given by 

G(x, = uj + /i — — x; (5) 

this implies a spectral function A(x,lo) = — -^ImG(x, uj), 
to be used later in the optical conductivity. 

Because of the simplicity of the non-interacting density 
of states, Eq. ([3]) can be integrated analytically. The 
result is: 



G(w) 



+ fi — E(w) 
2t 



2t 



1 



(6) 



where the y () is taken with a sign equal to sgn[Re(a; + 

For the SIAM, the impurity Green's function can be 
written as: 

1 



Gimp (w) 



+ fi-e d - A(w) - S imp (w) 



(7) 



in which A(u>) = X]fc=2 ^* /( w + M ~ e fc) is the hybridiza- 
tion function and Si mp is the impurity self-energy. The 
self-consistent process is based on the local nature of the 
quantum system in the limit of infinite dimensions, in 
which the on-site Green's function for the lattice model 
can be averaged over all momenta and only depends on 
the frequency as we obtain in Eq. Instead of di- 

rectly solving the dynamic Hubbard lattice model, we 
only need to solve the SIAM (with fewer degrees of free- 
dom) ; from this we obtain the impurity Green's function 
which should be the same as the on-site Green's function 
for the lattice: 



G(u>) = G imp (uj) 
E(w) = E imp (uj) 



(8) 
(9) 



Eqs. (|8I9I) are the self-consistency conditions. 

A solution is required for the impurity Green's func- 
tion. While an exact solution is available through QMG 
techniques, 8 Potthoff suggested a much faster though ap- 
proximate procedure known as the two-site DMFT, 10 
which he benchmarked for the Hubbard model. In this 
approach^ the SIAM with only two sites is diagonal- 
ized; one site is for the impurity and one site represents 
the bath, so n s = 2. Therefore the bath parameters are 
single numbers, as 6k — E c and Vk = V . Besides making 
the procedure significantly faster, the two-site DMFT is 
more transparent, as the self-consistency conditions arc 
analytic. The self-consistency conditions (|8I9I) are re- 
placed by conditions at high and low frequency, to give 
two new self-consistentcy conditions which relate directly 
to the bath parameters: 

riimp = n (10) 
V 2 = zA4 0) (11) 



where ni mp (n) is the filling for the impurity site (con- 
duction band in the lattice model). The parameter 
z = 1/(1 — Re(dS(0))/dcj) has the meaning of the quasi- 

particle weight in the metal phase and —J2<i j> tfj 
is the second moment of the non-interacting density of 
states. For the model adopted here, Eq. fl3J, this be- 
comes = t* 2 = 1. Therefore the right-hand-side 
of Eq. (fTTj) reduces to the quasiparticle weight. In fact, 
the procedure follows closely that given in Ref. [1(1 and 
the reader is referred to that publication for full details. 
Once the impurity problem is solved, one can obtain the 
density of states for the original lattice through 

A(uj) = — ImGfw + i0+) = p(u + n - ReCS(w))) (12) 

where the second equality only follows because the self 
energy is given by a two-pole approximation^ the self- 
energy is obtained from the self-consistent condition, Eq. 
© and Dyson's equation, 

E irop (cj + z0+) = Gg,(w + i0+f 1 - G imp (u> + t0+y\ 

(13) 

We will also be interested in the behavior of the optical 
conductivity, which within the same local approximation 
is given by the expression 8 

e¥V f +QO J /(e) -/(e + w) 

■KflfV J_ 00 CO 

+ 00 

dyp(y)A(y,e)A(y,e + ij), (14) 



where a is the lattice constant and v = a d is the volume 
of the unit cell in d dimensions. As stated above, the 
single particle spectral function, A(x,uj), is defined as 



A(x,u) = -~Im( — ■ \— 



(15) 



This function is immediately known once the self-energy 
is determined. 



III. RESULTS AND DISCUSSION 

A. Site Hamiltonian 

The Dynamic Hubbard model Hamiltonian, Eq. ([1}, 
consists of electron degrees of freedom that can move 
throughout the lattice, along with pseudospin degrees of 
freedom that reside at each lattice site. The pseudospins 
model the ability of the ions to 'react' to the different 
electronic configurations by changing the orbitals when 
electrons are and are not present. Following Ref. 0, we 
focus on the on-site Hamiltonian for electrons: 



H 



DHM 



uj Q a l x + guo<r\ + [U - 2guj al}m t n il , (16) 
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which is easily solved, given the presence of 0, 1, or 2 
electrons. Using the spin- 1/2 a z eigenstates, |+ >,|— >, 
as a basis, we find, that with n electrons present, the 
eigenstates are 



| n > a = u(n)\+ > +v(n)\- > 
\n> b = v(n)\+ > +u(n)\- > 



(17) 
(18) 



with eigenvalues: 



*~( n )a = 5 n , 2 U - u! Q y/l + g 2 



e{n) b = 5 nj2 U + cj ov /1 + g 2 (19) 
The eigenvector components are given by 



u(0) = u(l) = v(2) 
v(0) = «(l) = u(2), 



with 



1 



«"(o) = 

v 2 (0) = l(l + 



(20) 
(21) 



(22) 
(23) 



The expectation value of the pseudospin, in the ground 
state, illustrates the relaxation of the orbital required, 
depending on the number of electrons present. For ex- 
ample, the expectation value of a z , in the ground state, 
is given by 



< 0|a z |0 >=< l|<r,|l >= u 2 (0) - v 2 (0) = -g/^Jl + g 2 , 

(24) 

and 



< 2\a z \2 >= u 2 {2) - v 2 {2) = +g/y/l + g 2 . (25) 
Similarly, for the ground state, we obtain 



< 0|<r B |0 >=< 11(7,11 >=< 2\<j x \2 >= -l/Vl+5 2 . 

. (26) 

For large g the z-component of the pseudospin switches 
from close to —1 to a value close to +1 as the occupancy 
changes from one to two electrons, but does not change 
when the occupancy changes from zero to one electron. 
In contrast, the ^-component remains constant as the oc- 
cupancy changes. Note that this occurs independently of 
the value of the on-site Coulomb repulsion, U ; in par- 
ticular, the excitation energy associated with an excited 
pseudospin state is given by the difference of Eqs. (TlT))) , 



ft = 2w v 7 ! + g 2 - 



B. Mott transition 



(27) 



In this subsection, we examine how the dynamics of 
the auxiliary boson field affects the Mott transition in 
the Dynamic Hubbard model. In the next subsection we 
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FIG. 1: Bath parameters (a) E c and (b) V as a function 
of electron filling, for parameters with and without a Mott 
transition. For the static Hubbard model we use U = 8 and 
U = 4, as cases with and without a Mott gap, to illustrate 
the expected features in either case. The disappearance of 
V 2 at half-filling demonstrates the vanishing of quasiparticle 
spectral weight (see Eq. HI}), and the vertical increase of 
E c at half-filling is also a symptom of the Mott gap. Note, 
in the cases with non-zero g (dynamic Hubbard model), the 
asymmetry with respect to half-filling. 



consider effects related to electron-hole asymmetry and 
'undressing'. 

The effect of the dynamic field on the Mott transition 
is best addressed near half-filling. One of the signatures 
of this transition, as approached from the Fermi Liquid 
side, is the disappearance of quasi-particle weight (QW). 
As noted in Eq. (fTTj) and below, the self-consistent pa- 
rameter V measures the quasi-particle spectral weight. 
In Fig. 1 we show results for both self-consistent param- 
eters, E c and V for a number of parameters. Focusing 
on the static Hubbard model, it is clear that, within the 
DMFT approximation, a Mott transition takes place for 
a critical value of 4 < U c < 8, since, for U = 4 the 
parameter V remains non-zero over all electron fillings, 
while, for U = 8 the parameter V dips to zero at half- 
filling. Similarly, upon examining E c vs. n (Fig. 1(a)), 
one sees a vertical jump at half-filling (present for U = 8 
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FIG. 2: (Color online) Electron filling as a function of chem- 
ical potential, for the cases indicated in Fig. 1. Note the 
plateaus near half-filling, indicative of the Mott transition. 



but absent for U = 4) as the characteristic signature of 
the Mott phase. The other two parameter choices illus- 
trate that, for sufficiently large pseudospin coupling g, 
the Mott transition occurs, even for modest values of the 
bare Hubbard interaction, U. Thus the Mott transition 
is induced for U = 4 with g — 3.8, for example. This fact 
is further reinforced in Fig. 1(a), and redrawn in Fig. 2 
using the relation between E c and /i that is implicit in 
Eq. (fTU|) , where now a plateau is present near half- filling 
for the parameter set U — 4 and g — 3.8, thus indicating 
the occurrence of a Mott transition. 

Further insight into the occurrence of the Mott insu- 
lator (and the inherent particle-hole asymmetry in this 
model described in the next subsection) can be gained by 
examining the behaviour of the expectation values of the 
pseudospin variables. We plot in Fig. 3 the expectation 
values (a x ) and (a z ) as a function of filling for two sets 
of parameters, one in which a metal-insulator transition 
does not take place (g = 2, shown in pink), and one in 
which it does (g = 3.8, shown in green). We should note 
that our values are in quantitative agreement with those 
of Ref. (0) (note that our values for (<j x ) are negative 
while their's are positive). Fig. 3 makes clear that below 
half-filling the pseudospin expectation values are fairly 
constant as a function of electron filling. An expecta- 
tion value of (a z ) = — 1 maintains an effective U that 
is U + 2gu>0i much higher than U itself, so that double 
occupancy is restricted (see Fig. 5 below). 

A simple way to view the metal-insulator transition is 
through a variational approach, where, on the one hand, 
one uses the product state in which the electrons form a 
Fermi sea, while the pseudospins remain in their ground 
state, a Z i = —1, for all sites, i. A simple calculation 
yields 



FIG. 3: (Color online) The pseudospin expectation values 
< a x > and < a z > as a function of electron filling, for the 
dynamic Hubbard model, with g = 3.8 (with a Mott gap 
transition), and g = 2 (without a Mott transition). Note 
that the pseudospin plays very little role below half-filling, 
but undergoes a big change above half-filling. 



where n is the electron concentration, and we have used a 
simplifying assumption that the electron density of states 
is a constant, g(e) — l/(8t), appropriate to a 2D tight- 
binding model with nearest neighbour hopping t. This 
is essentially a Hartree calculation. The competing state 
is an insulator, with one electron per site (at least up to 
half- filling) , no hopping, and a similar pseudospin state, 
a z i — — 1, for all sites, i. The energy per site for this state 
is —gujQ. Therefore, restricting ourselves to half- filling, a 
metal-insulator transition will occur for U beyond U cr it, 
where U cri t = 8t — IgujQ. The critical value of U is clearly 
lower as g increases; this is because while the pseudospins 
remain in their ground state the effective value of on-site 
Coulomb repulsion is increased by the presence of the 
coupling to the pseudospin. Therefore U itself can be 
smaller and, in combination with the effect of the pseu- 
dospin, still instigate a metal-insulator transition. Note 
that the values of (a z ), as given by Eq. Q240 for low filling 
are —0.89 and —0.97 for g = 2 and g = 3.8, respectively; 
these are close to —1, as used in the variational calcula- 
tion, and also agree very well with the numerical results 
shown in Fig. 3. This in turn leads to a more accurate 
estimate of the bare on-site repulsion, given by^ 



C/ba 



u 



2g 2 uj 



(29) 



(H) 



Fermi 



N 



-2in(2 - n) - guj + (U + 2guj )n 2 /4, (28) 



which results when the background degree of freedom is 
not allowed to relax. 

As the filling increases above half-filling, the expecta- 
tion values of the pseudospins change markedly. Elec- 
trons are no longer able to avoid double occupancy, so 
the pseudospin steadily changes from the |— > state 
to the |+ > state, to lower the on-site energy from 
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[/ max = U + 2guj to U m i n = U — 2gu}Q- As Fig. 3 indi- 
cates, there is essentially a linear increase of < a z > from 
. 9 ~ — 1 to H — . 9 ~ +1 as the band becomes 

completely full (n = 2), as indicated by Eqs. (j2"4")l and 
(j25|) . The relaxation of the pseudospin degree of freedom 
results in a lower quasiparticle weight, as Fig. 1(b) indi- 
cates. In Fig. 2, the curve representing g = 3.8 (shown 
in green) approaches full occupation at n s» 4. This point 
can be understood by the fact that [/bare is approximately 
zero when n = 2 (see Eq. (|2"5)> ). so that the chemical 
potential goes to the top of the bare band (2t*) plus the 

=(<7 2 -l) » 1.7**. 



energy shift due to the pseudospin — = 

■ — 1 vV+i 

As pointed out in Ref. IjlOQ . for the Hubbard model, 
the two-site DMFT approximation correctly produces 
three peaks in the single particle density of states cor- 
responding to the lower and upper Hubbard bands, and 
a quasi-particle resonance peak at the Fermi energy, with 
quasi-particle weight z. For the dynamic Hubbard model, 
we show in Fig. 4 results for the parameters [7 = 4, 
g = 3.8 and uj = 0.5, along with results for the Hubbard 
model with U = 8, for a number of different electron 
densities. Because of the pseudospin degree of freedom, 
the spectrum should contain at least four peaks, due to 
the appearance of more than two poles in the self energy; 
this is clearly the case in Fig. 4. Below half-filling (Fig. 
4(a)), the peak structure is very similar to that found 
in the Hubbard model, with U — 8. At half-filling (Fig. 
4(b), n = 1.0), the resonant peak at the Fermi level has 
all but disappeared, indicative of the Mott transition. 
Also shown are the results for much higher filling (Fig. 
4(c), n — 1.5), where the results are clearly not sym- 
metric with those at n = 0.5 (Fig. 4(a)), and certainly 
no longer similar to the results for the Hubbard model 
with U — 8. As found in Ref. (7), the peak near the 
Fermi energy is considerably sharper at n = 1.5 com- 
pared with n = 0.5, indicating that this model is less 
free electron-like at high electron filling compared to low 
electron filling (as seen in Fig. 1(b) as well). 

Fig. 5 shows the double occupancy as a function of 
electron filling. Again, very little double occupancy is 
present below half filling, as expected, though it is clear 
that the presence of the electron-pseudospin coupling g 
suppresses the double occupancy near half-filling, and 
therefore enhances the Mott-like features of the Hubbard 
model. Above half-filling, the double occupancy quickly 
rises; though it is a more subtle effect here, the presence 
of g enhances the double occupancy since, as the filling 
increases the effective Hubbard U is decreasing due to 
the relaxation of the pseudospin degrees of freedom. 

In Fig. 6 we show the total energy as a function of 
electron filling, for a number of parameters, as indicated 
in the figure. For an empty or completely full band the 
numerical results agree with those given by the analytical 
results obtained for the site Hamiltonian, Eq. p^|) : 



E(n = 0) = -Luo^l+g 2 
E{n = 2) = U-u ^/T+ 



(30) 
(31) 
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FIG. 4: (Color online) A comparison of the spectral function 
for the static model with U — 8 to the one for the dynamic 
model with U = 4, g = 3.8, and ljo = 0.5, for three different 
fillings. The two models behave very similarly below half- 
filling; at half-filling both undergo a Mott transition with the 
disappearance of the width of the resonance peak at zero fre- 
quency. Above half-filling the two models differ markedly; in 
particular, the peak width at the Fermi level is considerably 
smaller for the dynamic Hubbard model compared with the 
static Hubbard model. 



The results in Ref. (0) (their Fig. 14) are also in excel- 
lent agreement with the exact results given by the above 
equations. 
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FIG. 5: (Color online) The dependence of double occupancy 
vs. electron filling, for various parameters, for both the static 
and dynamic Hubbard models. For sufficiently large g in the 
dynamic Hubbard model, the double occupancy is driven to 
zero at half-filling. 



FIG. 7: (Color online) The quasiparticle spectral weight, z 
vs. electron filling, for a variety of parameter values. For 
the dynamic Hubbard model, increasing the pseudospin fre- 
quency uuo leads to a steady decrease in the spectral weight, 
particularly at high value of electron filling, thus making the 
model more asymmetric with respect to half filling. 




FIG. 6: (Color online) The total ground state energy vs. elec- 
tron filling for the same sets of parameters as in Fig. 5. The 
values at n — and n = 2 agree with the analytically derived 
values, as explained in the text. 



FIG. 8: (Color online) Same as Fig. 7, except that now g is 
increased to a value g — 2. Increasing g also increases the 
asymmetry with respect to half filling. 



C. Electron- hole asymmetry and undressing 
phenomenology 

It was proposed in Ref. (0) that the dynamic Hubbard 
model describes superconductivity driven by "undress- 
ing" : namely, that when the Fermi level is near the top 
of the band, pairing of hole carriers will lead to transfer 
of spectral weight from high to low frequencies and in 
particular in an increase of the quasiparticle weight and 
a decrease in the effective mass. These effects should ap- 
pear both in the single particle spectral function and in 



two particle spectral functions such as the optical con- 
ductivity. 

In the present paper we do not consider pairing corre- 
lation functions, and hence we cannot ascertain from our 
results whether or not the model describes superconduc- 
tivity. However, we can study properties of the model 
under hole doping. It is expected that the effects dis- 
cussed in the previous paragraph should also occur both 
in the normal and the superconducting state for an al- 
most filled electron band as a function of increased hole 
doping 5 . 

For the parameters considered in the previous subsec- 
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FIG. 9: (Color online) Same as Fig. 7, but with a smaller 
value of U (here U =1), and we have added a higher fre- 
quency result (wo = 20). With respect to Fig. 7, Mott physics 
is de-emphasized; instead, a decrease in spectral weight as 
electron filling increases arises primarily because of the role 
of the pseudospin degree of freedom. Note that the curves 
arise from data obtained through converging the parameter 
V through the DMFT iterative process. The symbols come 
from integrating the spectral function peak at the Fermi level 
(see Fig. 12 below), and the good agreement is clear. 



tion, it was found that doping the full band with holes 
led to a decrease rather than an increase in the quasipar- 
ticle weight, in contradiction to these expectations. Thus 
we also do not expect superconductivity driven by "un- 
dressing" for those parameters. However, we find that 
the expected "undressing" behavior does occur for lower 
values of the on-site repulsion U and / or larger values of 
the coupling g as well as for higher values of the boson 
frequency luq. 

Figure 7 shows the behavior of the quasiparticle weight 
z (recall V = z 2 ) for U = 2, g = 1 and various values 
of uio versus band filling. It can be seen that for wo > 2 
the quasiparticle weight indeed increases when the full 
band is doped with holes (i.e. the quasiparticle weight 
decreases with electron filling as n —> 2). Fig. 8 shows 
the same behavior occurring already for ujq > 1 when we 
increase the value of the coupling to g = 2 with U = 2. 
Similarly, the same behaviour occurring for ujq > 1 can 
be obtained by reducing the on-site repulsion to U = 1 
while keeping g = 1 (Figure 9). For these parameter 
ranges a Mott transition at half-filling does not occur. 

As shown in Fig. 10, for large values of the frequency 
wo, the behavior of the quasiparticle weight versus band 
filling n is described approximately by the expression 

z=[l + (S-l^] 2 (32) 

with 

S 2 = — ^ (33) 
1+3 



FIG. 10: (Color online) Same as Fig. 9, except that now we 
have increased wo to very large values; in this limit we should 
recover the correlated hopping model, and the result at low 
hole doping (near full electron occupation) should follow the 
analytical result as indicated.— This is clearly the case. 




FIG. 11: (Color online) To disentangle spectral weight reduc- 
tion due to the Hubbard U vs. the dynamic pseudospin effect 
(g), we show the quasiparticle spectral weight z vs. electron 
filling n for several values of U, including U — 0. All numeri- 
cal results shown are for a large value of pseudospin frequency, 
wo = 50. Agreement with the analytical result is excellent for 
(7 = 0. Discrepancies for other values of U shown come pri- 
marily from the quasiparticle reduction due to U, which is 
not accounted for in the analytical result. 



as expected.— Also, as Fig. 10 shows, for smaller values 
of wo the n— dependence of z is qualitatively similar but, 
as n — > 2, the magnitude is larger than that given by 
Eq. (|32p . This dependence of quasiparticle weight on 
the boson frequency is consistent with the behavior found 
in Ref. 6 using a generalized Lang-Firsov transformation 
within an Eliashberg treatment. It was also found in that 
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work that smaller boson frequency enhances the tendency 
to pairing. 

Figure 11 shows that the analytical form, Eq. (j3"2")l 
is indeed very accurate when U is not present. Re- 
sults shown for increasing values of U indicate additional 
decreasing of quasiparticle weight that occurs due to 
well documented 'Hubbard physics'. Nonetheless, for all 
cases, the overall decreasing trend as a function of elec- 
tron filling is clearly coming from 'quasiparticle dressing' 
due to the pseudospin degree of freedom (Eqs. (|32|) and 



In Figure 12 we show the behavior of the single particle 
spectral function for one particular electron filling, n — 
1.8, for a variety of 'momenta'. Each spectral function 
consists primarily of two peaks, separated by roughly the 
pseudospin excitation energy flo given by Eq. (f27|) . The 
weight of the quasiparticle peak at the Fermi level, shown 
in the 5th panel, corresponds to the residue z plotted for 
all fillings in Fig. 9 (middle curve). As shown there, this 
quasiparticle weight decreases (increases) with increasing 
electron (hole) concentration. In Fig. 13 we show the 
behaviour of the optical conductivity obtained using Eq. 
(Tl4l) . for several electron fillings. An exact calculation 
of the optical conductivity for a two-site models shows 
that the optical conductivity can generally be divided 
into four contributions, as we now describe. 

There will be a Drude part centered at lj = with 
a width 1/r, normally due to elastic impurity scatter- 
ing. In our calculations, this part appears as a ^-function 
at the origin, with artificially imposed broadening (see 
Fig. 13). This component involves transitions between 
the two coherent parts of the spectral functions (see Eq. 
(|T4))h For low electron fillings these transitions give rise 
to essentially the entire frequency dependent conductiv- 
ity. 

Three more components contribute as the electron fill- 
ing comes close to and exceeds half-filling. They have 
characteristic frequencies w^, and weights C,;, for i — 
L,M,H, corresponding to 'low', 'medium', and 'high'. 
These designations are relevant only when U is suffi- 
ciently small, as in Fig. 13. Then the 'low' frequency 
part is peaked near lo l sa U/2 + v / (U/2) 2 + 4t 2 S 2 , and 
involves transitions between the ground state and an 
excited state with a doubly occupied site (both states 
have pseudospins in their ground states). The quantity 
S = (1|2) gives an estimate of the polaronic effect of the 
pseudospin excitation required as an electron undergoes 
a hopJ^ 

The characteristic frequency for the 'medium' range is 
lom « S1q = 2loq\J\ + g 2 . This is the excitation energy 
for a pseudospin (see Eq. (|2"T)l ). and transitions between 
states differing by a pseudospin excitation are responsible 
for this part of the conductivity. These transitions will 
play a more significant role as the electron-pseudospin 
coupling strength, g, increases. 

Finally, the high frequency characteristic frequency is 
ljh ss 2f^o, and requires a transition in which two pseu- 
dospin excitations are created. This necessarily tends 
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FIG. 12: (Color online) The spectral function for various en- 
ergies, as shown, for n — 1.8. Here, we have used (7 = 1, 
g = 1, and ojq — 2 (see pink curve in Fig. 9). The spectral 
functions consist primarily of two peaks, separated by an en- 
ergy corresponding to the pseudospin excitation energy, Qo- 
The 2nd last frame shows the spectral function at an energy 
corresponding to the Fermi level, and the weight under the 
peak at a; = corresponds to the quasiparticle residue, z. 
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to happen only when the electron concentration is very 
high, i.e. in the hole region. The weights of these var- 
ious contributions are difficult to estimate in advance; 
they will depend on the hopping overlap integrals. Fig. 
13 shows the optical conductivity at various fillings for 
U = l, .g = l, and ujq = 2.0; then lol ~ 2.0, u)m ~ 5.7 and 
uiji ~ 11.3. For these same microscopic parameters, we 
also show the expectation values of the z and x compo- 
nents of the pseudospin degree of freedom as a function 
of electron density in Fig. 14. Their values are given 
analytically at n = and n = 2: (<r x ) = + g 2 

for both limits and (a z ) = —g/ y/l + g 2 (g/ y/T+ g 2 ) for 
n = (n = 2), and DMFT obviously gets these correctly. 
The steady increase of (<r z ) with increasing electron den- 
sity reflects the increased occurrence of the excited pseu- 
dospin state in the ground state as the lattice becomes 
more crowded. 

A more in-depth understanding of the optical conduc- 
tivity comes from examining the spectral weight, parti- 
tioned into the various contributions, as described above. 
For concreteness we define the Drude conductivity to in- 
clude contributions in the range, < u> < ujl. The other 
contributions are defined as follows: low frequency in the 
range lol < w < ^>m, midrange for wjj < ui < uih, and 
the high frequency range for uj > ujh- This is shown 
in Fig. 15. First, note that the total spectral weight is 
generally asymmetric as a function of electron concentra- 
tion. Note, moreover, that as the electron concentration 
approaches zero or full filling, the spectral weight ap- 
proaches zero, since there are no carriers in either case. 
The Drude weight, defined as the area under the near- 
zero frequency portion of the conductivity, is also asym- 
metric; this is shown by the black triangular points. An 
integration over the Drude portion of Eq. (TT4"|) shows 
that the Drude spectral weight satisfies the sum rule^: 



(TT> T (u))dcd 
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2dh 2 v 
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and in the limit of infinite dimension d — > oo: 
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(35) 



Here ui* is the renormalized plasma frequency. We can 
see from Eq. (|3"51 that the Drude weight depends on 
the quasiparticle weight z and the DOS at ep = fx — 
i?e(E(0)). The result from Eq. {35]) is also plotted in 
Fig. 15, and is even more asymmetric than the total 
spectral weight. In addition, other contributions are also 
plotted as a function of electron filling. 

Note that the total weight is Drude-like both for n —> 
and for n — > 2, i.e. the contributions from the finite 
frequency portions fall off more quickly. Nonetheless, it 
is clear that spectral weight at high frequency is most 
intense in the hole-like region (n > 1) as compared with 
the electron-like region (n < 1). This can be understood 
by the following argument. First, we focus on the (minor) 
role of transitions involving the Hubbard U. These are 



described by the so-called 'low' frequency contribution, 
and they peak at half-filling. This is because that filling 
corresponds to the situation where more fluctuations are 
liable to occur. At low fillings U hardly plays any role, 
while at high filling U, by virtue of playing the same role 
for almost all electrons, again hardly plays any role. 

With excitations involving the pseudospin degree of 
freedom the situation is a little more subtle. At low to 
intermediate fillings, the pseudospins play almost no role 
because single electron transitions do not require a relax- 
ation of the pseudospin degree of freedom. At high to in- 
termediate fillings, however, transitions generally involve 
the pseudospin degree of freedom, because the overlap 
between a singly occupied state and a doubly occupied 
state, S = (1|2), can be much less than unity (in con- 
trast, for this model, the corresponding overlap between 
an empty and singly occupied site, T = (0|1) remains 
unity). Thus, there is an electron- hole asymmetry, with 
the pseudospin physics playing a large role in the elec- 
tron density region 1 < n < 2. In fact, variation in 
the pseudospin expectation value, (<J X ), reflects this fact 
also. Note the resemblance between the electron density 
dependence of (a x ) in Fig. 14 (inverted) with the elec- 
tron density dependence shown by the spectral weight 
in the medium frequency range shown by the blue as- 
terisks in Fig. 15 (the parameter regimes used in Fig. 
3 depict this variation of (a x ) in the hole-doped region 
even more pronouncedly). Hence, optical spectral weight 
in the medium frequency region plays a more important 
role for electron densities on the hole- like side of the phase 
diagram. This is consistent with the explanation in Ref. 
(0), that since (a x ) measures the fluctuations of o z (in 
imaginary time), it should achieve a minimum (they ob- 
tained a maximum for reasons that are unclear) at the 
point where (a z ) is changing the most quickly (as a func- 
tion of electron density). This tends to occur midway 
through the hole-like side of the phase diagram, i.e. near 
n rs 1.5. 

Such transfer of spectral weight from high to low fre- 
quencies as function of increasing hole concentration is 
seen in high T c cuprates, both in the single particle spec- 
tral function (photoemission experiments)^! and in the 
two particle spectral function (optical conductivity )iSr—. 



IV. CONCLUSIONS 

We have investigated various properties of the dynamic 
Hubbard model by using dynamical mean field theory 
in the two-site approximation. Where comparison was 
possible our results agree surprisingly well with Quantum 
Monte Carlo results 7 performed on small clusters (6 x 
6) at finite temperature (/3 = 5). In agreement with 
Ref. (0) we find that the presence of an auxiliary degree 
of freedom enhances the Mott transition at half filling. 
We see this in the vanishing of the quasiparticle spectral 
weight, in the appearance of a plateau in the n vs. /i 
curve, and in the vanishing width of a quasiparticle peak 
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in the density of states. 

An important property of this model is its asymme- 
try with respect to half filling. The hole side (n > 1) 
is always considerably more dressed than the electron 
side, because of the relaxation of the pseudospin degree 
of freedom. This occurs because electrons minimize their 
Coulomb repulsion, at a cost of becoming somewhat more 
sluggish in their movements, i.e. they form polaronic-like 
states. 

We have identified a parameter regime where the 
electron- hole asymmetry is very clear (see Figs. 7-11), 
and the quasiparticle spectral weight increases linearly 
with hole doping away from n = 2, as described in an 
effective model with correlated hopping, 2 This quasipar- 
ticle undressing is a general phenomenon that occurs not 
only as a function of hole doping, but as a function of, for 
example, temperature changes and phase transitions. 22 



Understanding the degree to which this undressing is ro- 
bust as the auxiliary degree of freedom is reduced from 
the anti-adiabatic limit to a more physical regime is one 
of the goals of this paper. Fig. 9 in particular illustrates 
that the degree of dressing (at n = 2) is reduced as uiq is 
reduced, and therefore the degree of 'undressing' as the 
electron occupation is decreased from n — 2 is reduced. 
Future work will determine the impact of this frequency 
scale on superconductivity. 

We computed single particle spectral weights and the 
optical conductivity for this model. In both cases the 
frequency dependence is distinctly different for electron- 
like and hole-like doping levels. The asymmetric behav- 
ior for holes vs. electrons is clear as a function of dop- 
ing; determining this asymmetry for a given doping level 
can be established through photoemission or tunneling 
experiments.— 
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FIG. 13: (Color online) The optical conductivity is shown at 
different filling for U = l,g = l,Wo = 2. Though not so 
apparent in this figure, there are four primary components 
with characteristic frequencies, uj — (Drude), ujl ~ U/2 + 
^(U/2) 2 + 4t 2 S 2 « 2, ui, = fl = 2CJ0V/H-3 2 « 5.7 and 
ujh = 2fi ~ 11.3. 
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FIG. 14: (Color online) The pseudospin expectation values 
< a x > and < <r z > as a function of electron filling with 
U = l,g = l,w = 2 
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FIG. 15: (Color online) The spectral weight of the opti- 
cal conductivity, as a function of electron filling. The total 
weight, along with the weight of each component (see figure 
caption in Fig. 13), is plotted as a function of electron fill- 
ing. As in Fig. 13, we have used the same parameter set 
U = l,g = l,u) = 2.0. 



